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ABSTRACT 

On galactic scales, the surface density of star formation appears to be well correlated 
with the surface density of molecular gas. This has lead many authors to suggest 
that there exists a causal relationship between the chemical state of the gas and its 
ability to form stars - in other words, the assumption that the gas must be molecular 
before star formation can occur. We test this hypothesis by modelling star formation 
within a dense cloud of gas with properties similar to a small molecular cloud using 
a series of different models of the chemistry, ranging from one in which the formation 
of molecules is not followed and the gas is assumed to remain atomic throughout, to 
one that tracks the formation of both H2 and CO. We find that presence of molecules 
in the gas has little effect on the ability of the gas to form stars: star formation can 
occur just as easily in atomic gas as in molecular gas. At low densities (< 10 4 cm -3 ), 
the gas is able to cool via C + fine-structure emission almost as efficiently as via CO 
rotational line emission, while at higher densities, the main cooling process involves 
the transfer of energy from gas to dust, meaning that the presence of molecules is 
again unimportant. Cooling by H2 is particularly inefficient, accounting for as little 
as 1 percent of the overall cooling in the cloud. Rather than the chemical makeup, 
we find that the most important factor controlling the rate of star formation is the 
ability of the gas to shield itself from the interstellar radiation field. As this is also a 
prerequisite for the survival of molecules within the gas, our results support a picture 
in which molecule formation and the formation of cold gas are both correlated with 
the column density of the cloud - and thus its ability to shield itself - rather than 
being directly correlated with each other. 

Key words: galaxies: ISM - ISM: clouds - ISM: molecules - stars: formation 



<3 



1 INTRODUCTION 

Essentially all Galactic star formation occurs within dense 
clouds of molecular gas, known as molecular clouds (MCs). 
Furthermore, recent observations have demonstrated that 
on large scales within local galaxies, there is a surprisingly 
close correlation between the surface density of star forma- 
tion, Esfr, and the surface density of molecular hydrogen, 
E h 2 (see e.g.|Wong fc Blitz 1 2002 1 |"Leroy et aI.|2008||Bigiel et 



al. 2008, 2011). An obvious interpretation of these observa- 



tions is that the formation of stars depends on the presence 
of molecular gas (|Krumholz & McKee|2005||Elmegreen|2007 



Krumholz et al. |2009|). However, it is not immediately ob- 



vious why this should be the case. Although H 2 can be an 
important coolant of interstellar gas (|Gnedin fc Kravtsov 



2011), it is effective only at temperatures of a few hundred 



Kelvin or higher. At the temperatures typical of gas within 
molecular clouds (T ~ 10-20 K), the H2 cooling rate is tiny, 
and hence H2 cannot play a direct role in enabling the gas 
to cool, undergo gravitational collapse, and form stars. 



A more promising route by which H2 can influence the 
thermodynamics of the gas is through the fact that its pres- 
ence is required for the efficient formation of carbon monox- 
ide (CO). Unlike H2, CO can provide effective cooling even 
in very low temperature gas, and it is generally found to 
be the dominant gas-phase coolant within prestellar cores 
( |Neufeld, Lepp fe Melnick|1995||Goldsmith|200l] . However, 
even if the H2 and CO were absent, the gas would still be 
able to cool to low temperatures through fine structure emis- 
sion from ionized and neutral atomic carbon. In previous 
modelling, we have shown that C + cooling alone can reduce 
the gas temperature to T ~ 20 K within gas that is shielded 
from the effects of photoelectric heating by dust extinctions 
of A v ~ 1-2 or more ( |Glover fe Mac Low||2007t . To re- 
duce the temperature further, to the 10 K typical of most 
prestellar cores, a molecular coolant such as CO is still re- 
quired, but it seems unlikely that the difference between the 
~ 10 K temperatures reachable with molecular cooling and 
the ~ 20 K temperatures reachable with atomic fine struc- 
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ture cooling can have any great effect on the ability of gas 
to form stars. 

One is therefore led to ask whether the presence of 
molecules is truly essential for star formation. Perhaps the 
reason that we find such a good correlation between molecu- 
lar gas and star formation is not that molecular gas is a pre- 
requisite for star formation, but instead that the formation 
of molecules and the formation of stars both correlate with 
some other factor, such as the existence of large, dense clouds 
of gas. In this picture, the molecules are a consequence, not 
a cause, of the conditions required to form stars. 

To help us better understand whether or not the pres- 
ence of molecular gas is a necessary condition for star forma- 
tion, we have performed a number of numerical simulations 
of the formation of stars in a dense cloud of gas with prop- 
erties similar to a small molecular cloud. We use a range 
of different models of the chemistry and thermal balance of 
the gas to explore the influence that molecule formation and 
dust shielding have on the ability of the cloud to form stars. 
In Section [2] of this paper, we discuss the numerical setup 
that we use to perform these simulations, and our choice 
of initial conditions. The results are presented in Section [3] 
and in Section [4] we discuss what these results imply for the 
nature of the link between molecular gas and star formation. 
We conclude in Section [5] 



2 SIMULATIONS 

2.1 Numerical approach 

Our simulations were performed using a modified version of 
the Gadget 2 smoothed particle hydrodynamics (SPH) code 
(Springel 2005). Our modifications include a sink particle al- 
gorithm for treating gravitationally collapsing regions that 
become to small to resolve (Bat e, Bonnell &: Price|1 995), the 
inclusion into the equation of motion of an optional confining 
pressure term ( |BenzfT990t , a treatment of gas-phase chem- 
istry (described in more detail below), as well as radiative 
heating and cooling from a number of atomic and molecular 
species ([Glover fe Jappsen||2007| |Glover et al.||2010| |Glover 



Table 1. Simulation details 



& Clark 201 1|, and an approximate treatment of the attenu- 



ation of the Galactic interstellar radiation field (ISRF) . The 
effects of magnetic fields are not included. 

In order to establish the role that molecular cooling 
plays in determining the star formation rate, we have run a 
series of simulations a small, dense gas cloud with proper- 
ties similar to those of a small molecular cloud, using several 
different treatments of the cooling and chemistry of the gas. 
In our least realistic model (run A), we include the simple 
network for hydrogen chemistry described in Glover & Mac 



Low 



(120071), together with the model for CO chemistry de- 
veloped by Nelson &; Langer (19971), but we assume that 
the gas remains optically thin throughout the simulation; 
i.e. we neglect any attenuation of the ISRF by dust absorp- 
tion or molecular self-shielding. In practice, this means that 
the gas remains in atomic form at all but the highest den- 
sities (see Section 



3.2 



below). In run B, we include the ef- 
fects of dust absorption, but ignore the chemical evolution 
of the gas, forcing it to remain in atomic form throughout. 
To model the dust absorption, we use our new Treecol algo- 
rithm, as described in more detail in Section [2.2| below. In 
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run C, we include the simple network for hydrogen chem- 
istry described in | Glover & Mac Low ( 2007 ) , but assume 

or 



that the carbon and oxygen in the gas remain as C 
O, respectively. In this model, the only molecular cooling 
comes from H2. Finally, in runs Dl, D2 and D3, we include 
bot h the[Gl over &; Mac Low] ( |2007| hydrogen chemistry and 
the Nelson & Langer (1997) CO chemistry, this time includ- 
ing the effects of dust absorption and self-shielding. This 
setup has been shown to produce more CO than would be 
the case with more detailed chemical treatments ( [Glover 



|Clark||201l| , and hence gives us a reasonable upper bound 
on the effectiveness of CO cooling. We performed three sim- 
ulations using this treatment of the chemistry: one in which 
the hydrogen is initially atomic (run Dl), and two in which 
it is initially molecular (runs D2, D3). These final two runs 
differ only in the numerical resolution adopted: run D2 is 
performed using our standard resolution (discussed in Sec- 
tion 2.4 below), while run D3 is performed using ten times 
more SPH particles, and hence has ten times better mass 
resolution. A brief summary of these models is given in Ta- 
ble CD 

None of the models that we study in this paper include 
the effects of freeze-out of species such as CO onto dust 
grains. However, we do not expect this simplification to sig- 
nificantly affect the thermal balance of the cloud (see e.g. 



Goldsmith||2001). 



2.2 Attenuation of the ISRF 

We assume that our model clouds are illuminated by the 
standard interstellar radiation field, using the parameteriza- 
tion from |Draine| p978| in the UV, and from |Black| ( [1994} 
at longer wavelengths. To model the effects of extinction, 
H2 self-shielding, CO self-shielding and the shielding of CO 
by H 2 , we use the Treecol algorithm (Clark et al., in prep.) 
This algorithm allows us to approximate the distribution of 
column densities (of hydrogen nuclei, molecular hydrogen 
or CO) seen by each particle by using information stored 
in the oct-tree that is also used in the Gadget 2 code to 
compute gravitational accelerations. Our current implemen- 
tation of Treecol uses the Healpix pixelation scheme ( |G6rski| 
et al.|2005 ) and yields for each SPH particle a 4tt steradian 
projection of the column density distribution seen by that 
particle, discretized into 48 equal- area Healpix pixels. 

To illustrate how we can use this pixelated projec- 
tion to compute photoheating or photodissociation rates 
for a given SPH particle, let us follow how we would go 
about computing the photoelectric heating rate. We be- 
gin with the appropriate pixelated column density distri- 
bution, which in this case is that for the column density of 
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Figure 1. Upper panel: mass in sinks as a function of time in runs 
B (shorted-dashed line), C (solid line), Dl (long-dashed line) and 
D2 (dot-dashed line). Middle panel: number of sinks formed in 
each simulation, plotted as a function of time. Lower panel: mean 
sink mass versus time in the same four simulations. For reference, 
the mass resolution in these simulations was 0.5 Mq. 



hydrogen nuclei in all forms, A^H,tot- We use the relation- 

-22/ at /-I nn% -2\ 



ship Ay = 5.348 x 10~^(jV H ,tot/l cm" 2 ) between jV H ,tot 
and the visual extinction Ay (Bohlin, Savage fe Drake|1978| 



Draine &; Bertoldi 1996|) to associate a visual extinction with 



each pixel. We next compute an attenuation factor / pe for 
each pixel, which we define as / pe = r p e(Av)/r p e(0), where 
r pe (iv) is the photoelectric heating rate in gas with a visual 
extinction equal to Ay, and r pe (0) is the photoelectric heat- 
ing rate in the optically thin limit. We take our expressions 
for these quantities from Bergin et al. (2004). Finally, once 
we have computed / pe for each pixel, we make use of the fact 
that the Healpix pixels have equal areas to compute a mean 
attenuation factor (/ pe ) for the gas represented by the SPH 
particle by averaging / pe over all pixels. The photoelectric 
heating rate in this gas is then just r pe = (/ pe )r pe (0). 

We use the same strategy to compute the rate at which 
dust is heated by the interstellar radiation field, taking our 
expression for the dependence of the heating rate on Ay 



from |Glover Sz Clark (2011 ). We also use a similar strategy 
to compute the H2 and CO photodissociation rates. How- 
ever, in this case it is necessary to use the Treecol algo- 
rithm to compute the distributions of H 2 and CO column 
densities in addition to the total hydrogen column density, 
since these are needed in order to compute the mean atten- 
uation factors corresponding to H 2 self-shielding, (/h 2 ,h 2 )j 
CO self-shielding, (/co.co), and the shielding of CO by H2, 
(/h 2 ,co)- Our expression for /h 2 ,h 2 comes from Draine & 
|Bertoldi (1996), while those for fu 2 ,co and /co,co come 
from Lee et al. (1996). A more accurate treatment of CO 



self-shielding and the shielding of CO by H2 has recently 
been given by Visser, van Dishoeck, &; Black (2009), but 
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Figure 2. Gas temperature plotted as a function of n, the number 
density of hydrogen nuclei, in runs B, C, Dl and D2 (panels 2-5) 
at a time immediately prior to the onset of star formation in each 
of these runs. Note that this means that each panel corresponds 
to a slightly different physical time. For comparison, we also plot 
the temperature of the gas as a function of density in run A (panel 
1, at the top) at a similar physical time, t = 2.3 Myr, although in 
this case, this is long before the cloud begins to form stars. 



we would not expect our results to change significantly if 



(1996) 



we were to use this in place of the older Lee et al. 
treatment. 



2.3 Sink particles 

Sink particles are formed once the gas exceeds a number 
density of 10 7 cm -3 , which is around the density at which 
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we start to lose resolution in the 2 x 10 6 particle simulations 
(assuming a minimum temperature of 8K; see Section 3.2 
below) . To ensure that sink particles are not created by mis- 
take, we require the parcel of gas comprising the candidate 
particle and its 50 nearest neighbours to be gravitationally 
bound and collapsing (following the prescription given by 
Bate, Bonnell & Price 1995). Once formed, sink particles 



3 RESULTS 



3.1 Star formation in the clouds 



are able to accrete additional SPH particles coming within 
a preset sink accretion radius r aC c, provided that the SPH 
particles are not only gravitationally bound to the sink, but 
are also more strongly bound to it than to any other sink 
within the computational volume. We prevent new sink par- 
ticles from forming within a distance of 2 r aC c of existing sink 
particles. In this study, we set r acc = 0.0026 pc, equal to the 
Jeans length at n — 10 7 cm -3 and T = 8 K when the mean 
molecular weight \i — 2.33 amu. All gravitational forces, 
including those from the sinks, are softened using the stan- 
dard kernel softening that has been shown by a number of 



authors to prevent 'artificial' fragmentation of the gas (Bate 



fe Burkert||1997| |Whitworth||1998| |Hubber et al.|2006| ). 



2.4 Initial conditions 

The starting point for our calculations is a uniform sphere 
with a hydrogen nuclei number density n = 300 cm -3 , and 
a total mass of 1O 4 M0, meaning that it has an initial radius 
of approximately 6 pc. Such conditions are typical of those 
found in nearby star- forming clouds, and give the cloud a vi- 
sual extinction at its centre of roughly three. We model the 
cloud with either 2 x 10 6 SPH particles (our standard reso- 
lution) or 2 x 10 7 SPH particles (our high resolution, used 
only in run D3), and we require that every SPH particle's 
quantities are calculated by averaging over 50 neighbours. 
This sets the mass resolutions at 0.5 and 0.05 Mq in our 



low and high resolution simulations respectively (Bate & 



Burkert 1997). We inject bulk (non-thermal) motions into 
the cloud by imposing a turbulent velocity field that has 
a power spectrum of P(k) oc k~ A . The energy in the tur- 
bulence is initially equal to the gravitational energy in the 
cloud, giving an initial root mean squared velocity of around 
3 kms -1 . The turbulence is left to freely decay via shocks 
and compression-triggered cooling. We set the initial tem- 
perature of the gas to 20K, but this is quickly altered when 
the simulation starts, as the gas tries to find an equilibrium 
between the various heating and cooling processes that are 
included in our model of the interstellar medium. We adopt 
a confining pressure equal to p e xt = 12000 K cm -3 , but note 
that runs with much smaller values of p ex t produce very 
similar results. 

We adopt total carbon and oxygen abundances relative 
to hydrogen that are xq — 1.4 x 10~ 4 and xq — 3.2 x 
10 -4 , respectively dSembach et dl. For the cosmic ray 

ionization rate of atomic hydrogen, we take a value £h = 
10 -17 s -1 . Cosmic ray ionization rates for other species (e.g. 
molecular hydrogen) were computed assuming that they had 
the same ratios to £h as given in the UMIST99 chemical 
database ( |Le Teuff, Millar fe Markwick||2000 ). 



In the top panel of Figure [T] we show the rate at which mass 
is converted into sink particles in runs B, C, Dl and D2. In 
all four of these runs, star formation begins within the cloud 
after roughly 2-2.5 Myr. For comparison, the global free- 
fall time of the cloud in its initial state is approximately 2.5 
Myr, and so in all four of these runs, stars begin to form 
within roughly one free-fall time. 

If we look in more detail at the time histories plotted 
in Figure^ we see that the time that elapses between the 
beginning of the simulations and the onset of star formation 
has a slight sensitivity to the chemical state of the gas. Run 
D2, which includes both H2 and CO chemistry, and which 
begins with the hydrogen already fully molecular, is the first 
run to begin forming stars. However, star formation in run 
C (initially atomic, H2 chemistry only) and run Dl (ini- 
tially atomic, H 2 and CO chemistry) is delayed only slightly 
compared to run D2, by roughly 0.1 Myr. Moreover, by the 
end of the simulation, all three simulations are forming stars 
at very nearly the same rate, although a greater number of 
stars have formed in run D2, owing to its head start. Finally, 
star formation in run B (atomic gas) is delayed for slightly 
longer, roughly 0.3 Myr, although once star formation does 
begin, it occurs at a very similar rate to that in the other 
three runs. 

Nevertheless, the differences between the star formation 
histories of the clouds simulated in these four runs are rela- 
tively small, despite the significant differences that exist in 
the chemical make-up of the clouds. The presence of H2 and 
CO within the gas appears to make only a small difference in 
the ability of the cloud to form stars. Furthermore, the fact 
that a cloud that does not form any molecules is not only 
able to form stars, but does so with only a short delay com- 
pared to one in which all of the hydrogen and a significant 
fraction of the carbon is molecular is persuasive evidence 
that the formation of molecules is not a prerequisite for the 
formation of stars. 

Examination of the nature of the stars formed in these 
four simulations is also informative. Naively, one might ex- 
pect that in the absence of molecular cooling, or more specif- 
ically CO cooling, the minimum temperature reached by the 
star-forming gas would be significantly higher. If so, then 
this would imply that the value of the Jeans mass in gas 
at this minimum temperature would also be significantly 
higher. It has been argued by a number of authors (e.g. 
| Larson|[2005l |Jappsen et"aT1[^505l |Bonnell, Clarke &: Bate] 
|2006 ) that it is the value of the Jeans mass at the mini- 
mum gas temperature in a star-forming cloud that deter- 
mines the characteristic mass in the resulting initial mass 
function (IMF). Following this line of argument, one might 
therefore expect the characteristic mass to be much higher 
in clouds without CO. However, our results suggest that this 
is not the case. In runs C, Dl and D2, the mean mass of the 
stars that form is roughly 1-1.5 M (Figure [I] bottom 
panel) and shows no clear dependence on the CO content of 
the gas. Indeed, the mean mass is slightly lower in run C, 
which has no CO, than in run D2, which does. The mean 
stellar mass that we obtain from these simulations is slightly 
higher than the characteristic mass in the observationally- 
determined IMF, which is typically found to be somewhat 
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less than a solar mass flChabrierpOOl) |Kroupa|2002) . How- 
ever, this is a consequence of our limited mass resolution, 
which prevents us from forming stars less massive than 0.5 
Mq, and hence biases our mean mass towards higher val- 



ues. (We return to this point in Section 3.3 below). In run 
B, we again find a greater difference in behaviour, but even 
in this case, the mean stellar mass remains relatively small, 
at roughly 2 M©. It therefore appears that the presence or 
absence of molecules does not strongly affect either the star 
formation rate of the clouds or the mass function of the stars 
that form within them. 

Conspicuous by its absence from our discussion so far 
has been run A, the run in which we assumed that the 
gas remained optically thin throughout the simulation. In 
this run, we find a very different outcome. Star formation is 
strongly suppressed, and the first star does not form until 
t — 7.9 Myr, or roughly three global free-fall times after the 
beginning of the simulation. The results of this run suggest 
that it is the ability of the cloud to shield itself from the 
effects of the interstellar radiation field, rather than the for- 
mation of molecules within the cloud, that plays the most 
important role in regulating star formation within the cloud. 



(c.f. Krumholz, Leroy fe McKee 2011) 



3.2 Thermal and chemical state of the gas 

3.2.1 Temperature distribution 

In order to understand why molecular gas appears to be of 
only very limited importance in determining the star for- 
mation rate, it is useful to look at the thermal state of the 
gas in the different runs at the point at which they begin 
forming stars. This is illustrated in Figure [2] for runs B, C, 
Dl and D2. For comparison, we also show the temperature 
distribution of the gas in run A at t = 2.3 Myr (i.e. at a 
similar time to the other four runs, albeit roughly 5.6 Myr 
before run A itself begins to form stars) . The first point to 
note is the basic similarity of the temperature distribution 
in most of the runs. In runs B, C, Dl and D2, the tempera- 
ture decreases from roughly 100 K at n ~ 10 cm -3 to 10 K 
at n = 10 5 cm -3 and 8 K at n = 10 6 cm -3 . This corre- 
sponds to a relationship between temperature and density 
that can be approximated as T oc p~ ' 25 at n < 10 5 cm -3 , 
or a relationship between pressure and density P oc p 0-75 , 
in good agreement with the relationship P oc p ' 73 proposed 
by [Larson] ( 1985, 2005). The fact that the effective equation 
of state of the gas is significantly softer than isothermal (i.e. 
P oc p) means that the local Jeans mass decreases rapidly 
with increasing density within the cloud, a factor which is 
known to greatly assist gravitational fragmentation (see e.g. 
Li, Klessen, fe Mac Low1[2003l |Jappsen et alT][2005| |Clark,| 



Glover fe Klessen|2008||Dopcke et al.|2011| ). In run A, on the 
other hand, there is a brief period of cooling at low densities, 
but the gas then begins to reheat at densities n > 500 cm -3 , 
and the temperature in the densest gas is higher than the 
temperature at the edge of the cloud. It is therefore not 
surprising that gravitational fragmentation is strongly sup- 
pressed in this run. 

If we look in more detail at the relationship between 
temperature and density in the four runs that do show sig- 
nificant cooling, then several additional features become ap- 
parent. First, all of the runs show a significant degree of scat- 



ter in the relationship at densities between n ~ 100 cm -3 
and n ~ 10 4 cm -3 (in run B) or n ~ 10 5 cm -3 (in the other 
runs), but this scatter abruptly vanishes at higher densities. 
This pronounced change in behaviour is due in part to a 
change in the dominant coolant within the clouds. At low 
densities, gas-phase coolants (primarily C + and CO) domi- 
nate, while at higher densities, energy transfer between gas 
and dust becomes the dominant process regulating the gas 
temperature. In the regime dominated by gas-phase cooling, 
the local cooling rate depends on the local value of the opti- 
cal depth in the appropriate cooling line, and hence on the 
local details of the velocity structure of the cloud. In the 
dust-dominated regime, on the other hand, the gas temper- 
ature quickly converges to the dust temperature, which is 
determined largely by continuum emission and absorption 
and which is thus insensitive to the velocity structure of the 
cloud. In addition, another effect contributing to the scatter 
at low densities is the sensitivity of the photoelectric heating 
rate to the visual extinction. We know from previous work 
that the mean visual extinction within a given fluid element 
in a turbulent cloud is only poorly correlated with the vol- 
ume density of the cloud ( Glover et al.|20To| , and hence the 
value of the photoelectric heating rate shows a significant 
scatter at any given density. At low densities, photoelectric 



heating is the dominant heat source (see Section 3.2.4 be- 
low), and so this scatter in the heating rate helps to create 
a significant scatter in the gas temperature. At higher den- 
sities, photoelectric heating becomes far less important, and 
so the scatter in the photoelectric heating rate has much less 
effect on the temperature. 

A second feature to note is the pronounced "hump" 
in the temperature-density distribution, roughly centered at 
T rsj 20 K and n~2x 10 4 cm -3 , that is present in runs C 
and Dl, but absent in runs B and D2. This feature is caused 
by H2 formation heating: we assume (following Takahashi & 
|Uehara|2001| that H2 molecules form on dust grains with a 
high degree of rotational and vibrational excitation, and at 
densities above a few thousand particles per cubic centime- 
tre, much of this energy is converted into thermal energy by 
collisional de-excitation of the newly- formed H2 molecules. 
This does not occur in run B because no H2 is allowed to 
form in that run, and it does not occur in run D2 because the 
hydrogen starts in fully molecular form, and although some 
is subsequently dissociated, the atomic hydrogen fraction at 
the appropriate densities remains small (see Figure [4]). 

Finally, the influence of the CO in runs Dl and D2 is 
clear if we look at the temperature distribution of the gas at 
densities of around 10 3 to 10 4 cm -3 . The minimum gas tem- 
perature reached in the runs with CO is approximately 5 K, 
roughly half the size of the minimum temperature reached 
in the runs that must rely on C + cooling. However, despite 
the fact that the gas is able to reach lower temperatures in 
the runs with CO cooling, very little of it does so: most of 
the gas at these densities has a temperature that is close 
to 10 K, implying that the difference in the mean temper- 
ature of the gas is much smaller than the difference in the 
minimum temperature. 

Given the similarity between the temperature distribu- 
tion in run B and that in runs C, Dl and D2, it may at first 
seem odd that star formation in the former is delayed by 
roughly 0.3 Myr compared to star formation in the latter, 
and that the mean mass of the stars that form is noticeably 
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larger. However, this has a simple explanation. At a given 
density and temperature, the Jeans mass scales with the 
mean molecular weight \i as Mj oc /i _3//2 . In run B, the gas 
is fully atomic at all densities, and so \i — 1.27 amu (where 
we have assumed a helium fraction of 0.1 by number, rela- 
tive to hydrogen). On the other hand, in runs C, Dl and D2, 
the hydrogen in the densest gas is almost fully molecular, 
as illustrated in Figure [4] and so \i — 2.33 amu. Therefore, 
for any given choice of density and temperature within the 
dense gas, the Jeans mass is roughly a factor of 2.5 times 
larger in run B than in runs C, Dl or D2. If we assume that 
gas is assembled into dense filaments and cores at roughly 
the same rate in all of these runs, which seems likely given 
the similar temperature-density relationship, then this im- 
plies that in run B it will take longer for these structures 
to grow to the point at which they become gravitationally 
unstable, and hence that the onset of star formation will 
be delayed. It also suggests that the stars that do form will 
have systematically larger masses, in agreement with what 
we find in our simulations. 



3.2.2 Density distribution 

We can gain additional insight into the physical processes 
determining the star formation rate in our model clouds by 
examining the probability density function (PDF) of gas 
density generated in the different runs. In Figure |3j we show 
the mass- weighted density PDF of the gas in runs B, C, Dl 
and D2 immediately prior to the onset of star formation in 
these four runs. For comparison, we also show the density 
PDF in run A at a similar physical time, t — 2.3 Myr. 

In the four star- forming runs, the density PDF is broad, 
and shows evidence for a power- law tail at high densities, 
n > 10 4 cm -3 . Similar features have been found in pre- 
vious simulations of self-gravitating supersonic turbulence 
(e.g. |Klessen||20QQ| |Federrath et a"T1|2008| [Kritsuk, Norman 



fe Wagner||2011| ), in which run-away gravitational collapse 
has occurred. There is also now observational evidence for 
the presence of power-law tails in the density PDFs of ac- 
tive star-forming regions in real molecular clouds ( |Kainu- 
lainen et al"1|2009| |Lombardi, Lada Alves||2010| |Kainu- 
lainen et al.||2011). In contrast, the density PDF in run A 



is much narrower, and appears to have a log-normal form, 
as would be expected in a cloud in which gravitational col- 
lapse has not yet occurred. The higher gas temperatures in 
run A compared to the other runs mean that the turbu- 
lence has a much smaller root-mean-squared Mach number, 
and hence produces weaker compressions, leading to a sig- 
nificantly narrower PDF (|Padoan, Nordlund &; Jones||1997| 
Federrath, Klessen, fe Schmidt||2008| |Padoan fe Nordlund] 
2011). As a result, none of the gas in run A has yet reached 
a high enough density to become locally self-gravitating and 
to decouple from the larger-scale flow, while in the other four 
runs, a significant fraction of the gas mass is found in col- 
lapsing, self-gravitating structures. 

The fact that the four star-forming runs have very simi- 
lar density PDFs is a consequence of the fact that they have 
very similar temperature distributions. It is therefore easy 
to understand why they form stars at very similar rates: in 
all four cases, the amount of dense gas available for star for- 
mation is very similar, and the rate at which this dense gas 
is converted into stars is limited primarily by the local free- 



0.01 




n [cm 3 ] 

Figure 3. Comparison of the mass-weighted density PDF in the 
various runs. For runs B, C, Dl and D2, we show the state of the 
gas immediately prior to the onset of star formation, while for 
run A, we show the state of the gas at a similar physical time, 
t = 2.3 Myr. 



fall time, which does not depend upon the nature of the gas 
coolants. 



3.2.3 Molecular abundances 

In Figure [4] we plot the fractional abundance of H2, xh 2 , 
as a function of the hydrogen nuclei number density n in 
runs C, Dl and D2 at a time shortly before the onset of star 
formation in these simulations. For comparison, we also plot 
xh 2 versus n at a similar output time in run A. We show no 
results from run B because xh 2 = by design in this run. 
We define xh 2 as the ratio of the H2 number density nu 2 and 
the hydrogen nuclei number density n, which means that in 
fully molecular gas, xh 2 = 0.5n/n = 0.5. 

If we compare the results from runs C, Dl and D2, we 
see that in each case, the cloud is predominantly molecular 
at the point at which it forms stars. The mean H2 frac- 
tional abundance falls significantly below 0.5 only for den- 
sities n < 100 cm -3 , corresponding to warm, unshielded gas 
at the boundaries of the clouds. The rapidity with which 
H2 forms within the clouds is a consequence of the many 
transient density enhancements produced by the supersonic 
turbulence, which boost the overall H2 formation rate in the 



cloud by a factor of a few (for more details, see e.g. Glover & 
Mac Low|[2007| |Milosavljevic et al.||20TTf ). The main differ- 
ence of note between runs C, Dl and D2 is the larger scat- 
ter in xh 2 at high densities in runs C and Dl than in run 
D2. This reflects the difference in initial conditions between 
these runs: runs C and Dl start with all of their hydrogen in 
atomic form, and although most of it has become molecular 
by t rsj 2.5 Myr there are still a few regions that retain a high 
atomic fraction. In run D2, on the other hand, the hydro- 
gen is initially in fully molecular form, and atomic hydrogen 
is produced only by photodissociation. The main conclusion 
that we can draw from this is that at the point at which star 
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Figure 4. Fractional abundance of H2 (defined as the ratio of 
the H2 number density wh 2 and the number density of hydrogen 
nuclei, n), plotted as a function of n, for runs A, C, Dl and D2. 
With this definition of the fractional abundance, a value of 0.5 
corresponds to fully molecular gas. For runs C, Dl and D2, we 
show the state of the gas immediately prior to the onset of star 
formation, while for run A, we show the state of the gas at a 
similar physical time, t = 2.3 Myr. 



formation begins, the clouds have not yet reached chemical 
equilibrium throughout their volume. 

In run A, on the other hand, the situation is very differ- 
ent. In the absence of H2 self-shielding and dust shielding, 
the H2 fractional abundance is much smaller: the peak value 
at the displayed output time of 2.3 Myr is xh 2 — 2 x 10 -3 , 
implying that even in the densest gas, less than 1% of the 
hydrogen is in molecular form. In addition, the very small 
degree of scatter in the relationship between xh 2 and n im- 
plies that the gas is in chemical equilibrium, as we would 
expect given that in the absence of shielding, the H2 pho- 
todissociation timescale t p d — 600 yr, which is very much 
shorter than the age of the cloud. 

In Figure [5] we show how the fractional abundance of 
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Figure 5. As Figure^] but for the fractional abundance of CO. 
In this case, results are plotted only for runs A, Dl and D2, as 
the CO content of the gas in runs B and C is zero, by design. 



CO, xco = rico/n, varies as a function of n in runs A, Dl 
and D2. Results are shown at the same output times as in 
Figure^] No results are plotted for runs B or C, because 
xco = by design in these runs. Several points should be 
noted. First, the CO abundance in run A is tiny, as is to 
be expected given the low H2 abundance and the absence of 
dust shielding. Second, the distribution of CO abundances 
in runs Dl and D2 are very similar, despite the difference 
in the initial conditions for these two runs. This suggests 
that the CO content of the gas is primarily sensitive to its 
current H2 content and density structure. Any sensitivity to 
the chemical history of the gas, which is very different in runs 
Dl and D2, must be small. Additional support for this point 
is provided by Figure [6] which shows a CO column density 
projection of the clouds in runs Dl and D2 at t — 2.2 Myr. 
Although there appears to be somewhat more structure at 
low CO column densities in run D2 than in run Dl, the main 
features of the cloud are very similar in both simulations, 
and there is no clear observational signature of the different 
chemical histories that one could point to in these runs. 

3.2.4 Dominant thermal processes 

To help us to understand why the presence or absence of 
molecules appears to have such a limited effect on the tem- 
perature structure of the clouds and the progress of star 
formation within them, it is useful to examine how the dom- 
inant heating and cooling process vary with density within 
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Figure 6. CO column density projection at t = 2.2 Myr in run Dl (left-hand panel) and run D2 (right-hand panel). More small-scale, 
low column density structure is visible in run D2, which started with its hydrogen in molecular form, than in run Dl, which started with 
its hydrogen in atomic form. However, the high column density structures are quite similar in both cases. 



an atomic cloud and a molecular cloud. We have therefore 
calculated the median heating/cooling rate as a function of 
density for each process included in our thermal model for 
runs B and D2 at a time just before the onset of star for- 
mation. The results are plotted in Figure [7] (for run B) and 
Figure [8] (for run D2). Values are plotted only for the most 
important processes; minor contributions to the cooling from 
e.g. H + recombination or the collisional ionization of atomic 
hydrogen are omitted. We have chosen to plot the median 
values rather than the mean values because the latter are 
strongly influenced by dissipation in shocks, which occupy 
a relatively small fraction of the simulation volume but dis- 
sipate a significant amount of energy within the shocked 



sion, but the presences of H2 also introduces additional heat- 



ing processes such as H2 photodissocation heating ( | Black fe 
Dalg arno|1977 ), or heating due to the pumping of highly ex- 



regions (see e.g. |Pan &; Pad oan 2009). The median values 
give a fairer picture of the balance of heating and cooling 
rates in the bulk of the gas. 

Figure [7] demonstrates that in the atomic cloud, there 
are three main regimes. At densities n < 2000 cm -3 , the 
heating of the gas is dominated by photoelectric emission 
from dust grains, and the cooling is dominated by fine 
structure emission from C + . Between n ~ 2000 cm -3 and 
n~2x 10 4 cm -3 , C + remains the main coolant, but pho- 
toelectric heating becomes much less effective, owing to the 
increasing visual extinction of the cloud at these densities. 
In this density range, adiabatic compression of the gas (indi- 
cated in the plot as the pdV term), dissipation of turbulent 
kinetic energy in shocks, and cosmic ray ionization heating 
are all significant sources of heat. Finally, at densities greater 
than n~2x 10 4 cm -3 , the timescale for energy transfer be- 
tween gas and dust becomes short enough to couple the gas 
and dust temperatures together, at which point this (fol- 
lowed by thermal emission from the dust grains) becomes 
the most important cooling process. Weak shocks and adia- 
batic compressions together dominate the heating of the gas 
in this regime, each contributing close to half of the total 
heating rate. 

Figure [8] shows that in the molecular cloud, a number of 
additional processes come into play: H2 and CO provide ad- 
ditional cooling channels through their rotational line emis- 



cited vibrational levels of H 2 by UV photons (Burton, Hoi- 
|lenbach &; Tielens 1990 ) . Nevertheless, despite the additional 



complexity, we can again identify three main regimes marked 
out by different dominant processes. At n < 1000 cm -3 , the 
behaviour is very similar to that in the atomic run: C + is 
the dominant coolant, while most of the heating comes from 
photoelectric emission from dust. Above n = 1000cm -3 , C + 
quickly gives way to CO as the dominant coolant, reflecting 
the fact that the gas becomes CO-dominated at around this 
density (see Figure [5J , and the photoelectric heating rate 
also begins to fall off with increasing density. The fact that 
these changes occur at a very similar density is no coinci- 
dence: the photoelectric heating rate and the CO photodis- 
sociation rate have a very similar dependence on the visual 
extinction of the gas, and so both become unimportant at 
roughly the same point. CO remains the dominant coolant 
between n = 1000cm -3 and n ~ 10 5 cm -3 , but photoelectric 
heating quickly becomes irrelevant, and dissipation in shocks 
becomes the main source of heat. Finally, at n > 10 5 cm -3 , 
dust takes over from CO as the most important coolant, and 
pdV heating becomes almost as important as shock heating. 

Figure [8] also illustrates that cooling by H 2 is never par- 
ticularly important: at best, it contributes only a few percent 
of the total cooling rate, and at most densities contributes 
far less than this. In addition, it demonstrates that H 2 for- 
mation heating is unimportant in run D2, which is unsur- 
prising given the fully molecular initial conditions used for 
this run. A similar plot for run Dl would show a much larger 
contribution from H 2 formation heating at densities between 
n = 10 3 cm -3 and n - 10 4 cm -3 . 

If we compare Figures [7] and [8] we can see why the 
presence of H 2 and CO appears to have such a limited effect 
on the behaviour of the cloud. Below n — 1000 cm -3 (which 
is, let us not forget, more than three times higher than the 
initial mean density of the gas) , the cooling is dominated by 
C + in both simulations, and so H 2 and CO have very little 
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influence on the thermal behaviour of the gas. At higher 
densities, CO takes over from C + as the dominant coolant 
in run D2. However, if we compare the CO cooling rate at 
e.g. n — 10 4 cm -3 in run D2 with the C + cooling rate at 
the same density in run B, we see that they are surprisingly 
similar - the C + cooling rate is smaller than the CO cooling 
rate, but only by a factor of 2-3, despite the much larger 
energy required to excite the C + line compared to the CO 
J — 1 — >• line. It is this fact that ultimately renders the 
CO inessential: if we were to remove it, the temperature of 
the gas would rise slightly, and the cooling rate would fall by 
a factor of a few, but otherwise, the behaviour of the cloud 
would hardly change, and hence its ability to form stars 
would not be greatly affected. Since, as we have already seen, 
H2 cooling is also unimportant, we can therefore conclude 
that molecular cooling in general is not required in order to 
form stars, given the initial conditions examined here. In a 
realistic model, the gas that forms stars will be molecular 
(as in runs C, Dl and D2), but even if it were not (as in run 
B), it would still form stars at almost the same rate. 

Finally, we note that although the chemical model used 
in our simulations does not include atomic carbon, which 
could potentially be an important low-temperature coolant, 
we would not expect its inclusion to significantly affect the 
star- for mat ion rate. Atomic carbon is unlikely to provide 
more cooling at low temperatures than CO, and as we have 
already seen, the star formation rate is insensitive to the 
amount of CO present in the gas, implying that the star 
formation rate will also be insensitive to the ratio of C + 
to C. We have tested this assumption by resimulating the 
cloud from run D2 using a more detailed chemical model 
that does include atomic carbon (the NL99 model described 
in |Glover &; Cl ark 2011] which is a slightly modified version 
of a model first introduced by Nelson & Langer 1999). As 
expected, we find that the star formation rate in this run is 
essentially the same as in run D2. 



3.3 Influence of numerical resolution 

To establish the extent to which our results are affected 
by the limited mass resolution of our simulations, we per- 
formed a single simulation with the same initial conditions 
and chemical model as run D2, but which had ten times more 
SPH particles, and hence a mass resolution of 0.05 M©, ten 
times smaller than in our standard runs. We denote this run 
as D3. In Figure [5] we compare the growth of the mass in- 
corporated into sink particles, the number of sink particles 
formed, and the mean sink particle mass in runs D2 and D3. 
We see that the high resolution run D3 begins to form stars 
at almost the same point as the lower resolution run D2, and 
that the stellar mass increases at very nearly the same rate, 
although the much higher computational cost of the higher 
resolution run means that we are unable to follow the evo- 
lution of the system for as long as in the low resolution run. 
We also note the high resolution run forms a significantly 
larger number of stars than the low resolution run, with a 
mean stellar mass that is a factor of a few smaller. 

These results suggest that the rate at which stars form 
within our model clouds is not strongly influenced by the 
limited mass resolution of our simulations. This is what one 
would expect if the main factor determining the star forma- 
tion rate is the formation of self- gravitating, dense cores (see 
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Figure 9. Upper panel: mass in sinks as a function of time in runs 
D2 (dot-dashed line) and D3 (dashed line). These two runs use 
the same chemical model and initial conditions, and differ only in 
numerical resolution, with run D3 having ten times better mass 
resolution than run D2. Middle panel: number of sinks formed in 
each simulation, plotted as a function of time. Lower panel: mean 
sink mass versus time in the same four simulations. For reference, 
the mass resolution in run D2 was 0.5 Mq, while in run D3 it was 
0.05 M . 



e.g. |Krumholz fe McKee|[2005l |Padoan fe Nordlund|[20TT| , 
as most of these cores are well resolved even in our low res- 
olution simulations. The mass distribution of the stars that 
form, on the other hand, is sensitive to fragmentation occur- 
ring within the cores, which is not well resolved within our 
low resolution models, as the differences in the lower two 
panels of Figure [9] clearly demonstrate. 

We therefore conclude that our major results, which 
concern the timing and rate of star formation within the 
different clouds, should be insensitive to our choice of nu- 
merical resolution, but that higher resolution simulations 
will be needed to fully explore the sensitivity of the stellar 
initial mass function to the chemical make-up of the gas. 



4 DISCUSSION 

Our simulations demonstrate that once a dense, gravitation- 
ally unstable cloud of gas has formed in the ISM, molecu- 
lar cooling is not required in order to convert this gas into 
stars. In the absence of molecules, the cooling provided by 
the C + fine structure line at n < 2 x 10 4 cm -3 and by dust 
at n > 2 x 10 4 cm -3 are sufficient to allow stars to form. 
However, we cannot conclude solely from this that molec- 
ular gas is an inessential component of the star formation 
process. A skeptic could well argue that molecular cooling 
may be required in order to form the dense cloud of gas that 
we take as the initial conditions for our simulations, an as- 
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sertion that our current results can obviously do nothing to 
address. Fortunately, the issue of the assembly of cold, dense 
clouds of gas from the warm, neutral ISM has already been 
investigated by a number of authors (see e.g. |Koyama fc In-| 
utsuka 2002; Audit & Hennebelle 2005; Vazquez- Semadeni 
et al. H20071 IHeitsch et al.| |20081 |Heitsch fe Hartmann||2008T 



Banerjee et al.||2009 ) . 



These studies typically make use of a setup that in- 
volves converging flows of warm atomic gas, with an initial 
temperature and density corresponding to the warm neu- 
tral medium (WNM). The thermal models adopted vary 
somewhat in complexity, from a simple two-parameter cool- 



ing function introduced by Koyama & Inutsuka (2002) and 



used in a number of later works (e.g. Vazquez- Semadeni et 



al. 



2007 



Banerjee et al. 2009|, to the more detailed treat- 



ment used in e.g. Audit & Hennebelle (2005) that accounts 



explicitly for photoelectric heating from dust and fine struc- 
ture cooling from C + and O. However, one common fea- 
ture is that these models do not account for cooling coming 



from H2 or CO. Indeed, they typically do not include any 
non-equilibrium chemistry. Despite this, these models read- 
ily produce cold, dense clouds of gas in the region where the 
converging flows intersect. These cold clouds are produced 
initially by a thermal instability acting in the WNM that 
relies primarily on C + cooling. The clouds can then collide 
and agglomerate, forming gravitationally bound structures. 

The typical densities and masses of the cold clouds 
formed in these converging flow studies span a wide range of 
values, depending on the size of the simulation volume and 
the effective resolution of the simulation. Large-scale con- 
verging flow simulations, such as that performed by |Baner^] 
jee et al. (2009), have shown themselves capable of produc- 



ing cold clouds with properties similar to those we take as 
our initial conditions (i.e. a mean density of a few hundred 
cm -3 and a cold gas mass of 1O 4 M ) even though these sim- 
ulations do not include molecular cooling. Models of cloud 
formation on the scale of the Galactic disk also appear to 
have no difficulty in producing large cold gas clouds with- 
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out needing to invoke molecular cooling (see e.g. Tasker & 
Tan 2009). Taken together with our results, this implies that 
molecular gas is not necessary for star formation: it is pos- 
sible to form both cold, dense clouds, and also stars within 
these clouds, without needing to rely on molecular cooling. 

If molecular gas is not a prerequisite for star formation, 
then why is it that we find such a good observational cor- 
relation between the surface density of molecular gas, and 
the star formation rate surface density? An obvious expla- 
nation is that both the presence of molecular gas and the 
propensity of a gas cloud to form stars are correlated with 
some third factor, such as the column density of the cloud. 



An argument along these lines was given by Schaye ( 2004 ) , 



who pointed out that in regions of the Galactic disk with 
gas surface densities of less than 10 M pc -2 , ionization and 
heating of the gas by the extragalactic radiation field pre- 
vent the formation of a cold, neutral phase. The transition 
from the warm neutral medium to the cold neutral medium 
is associated with a dramatic decrease in the Jeans mass of 



the gas, and Schaye (2004) argues that this is the trigger 



for star formation. The transition from atomic to molecular 
gas coincidentally occurs at roughly the same surface den- 
sity (see e.g. Browning, Tumlinson & Shull 2003; Krumholz 
et al.|[2009 ), and so the result is a correlation between the 
presence of molecular gas and the formation of stars. 



The Schaye (2004) model assumes that all cold clouds 



form stars, but as we have seen in this paper, this is proba- 
bly an oversimplification. In a more recent paper, [Krum holz , 
Leroy fe McKee| ( |2011| ) look in more detail at the chemistry 
and thermodynamics of cold clouds in the ISM. They use 
simple ID cloud models that assume chemical and thermal 
equilibrium to explore a range of different cloud densities 
and visual extinctions, and show that in these models, the 
transition from atomic to molecular hydrogen is well cor- 
related with a further sharp drop in the equilibrium gas 
temperature. This correlation is not a result of H2 cooling. 
Rather, it occurs because the conditions required in order 
to attain a low gas temperature - high densities to boost 
C + cooling, and efficient dust shielding to suppress photo- 
electric heating - are similar to those required to produce a 
high equilibrium H2 fraction. They then argue that star for- 
mation is strongly correlated with regions of cold gas, owing 
to the T 3 / 2 temperature dependence of the Bonnor-Ebert 
mass scale, which makes gravitational fragmentation much 
easier to bring about in cold gas than in warm gas. In both 
of these models, the correlation between H 2 and star for- 
mation comes about because the H2 traces (but does not 
cause) the regions where the thermal pressure is low enough 
to allow stars to form. 

The results from our present study provide strong sup- 
port for this picture. They show that H2 cooling plays an 
insignificant role in determining the cloud temperature, that 
the differences between the temperatures of clouds cooled 
solely by C + or by a mix of C + and CO are very similar, 
and that the key factor enabling star formation within the 
clouds is the shielding of the interstellar radiation field by 
dust. In the absence of this shielding, the temperature of the 
gas remains high, star formation is strongly suppressed, and 
the H 2 abundance is very small. 



5 CONCLUSIONS 

In this study, we have investigated whether or not the forma- 
tion of molecular gas is a prerequisite for star formation. We 
have performed simulations using several different chemical 
models: one in which the gas is assumed to remain atomic 
throughout, a second in which H2 formation is included, but 
CO formation is not, and a third which follows both H2 and 
CO formation. We find only minor differences in the nature 
and rate of star formation in these simulations. In contrast, 
disabling the effects of dust shielding has a very strong ef- 
fect: the gas temperature does not fall much below 100 K 
and the formation of stars is strongly suppressed. 

We infer from these results that the observational cor- 
relation between H2 and star formation is not a causal re- 
lationship: H2 and CO are not required for star formation, 
and the fact that we find a good correlation between the 
H2 surface density and the star formation rate surface den- 
sity simply reflects the fact that both are correlated with 
some third factor. Our results suggest that the key factor 
is the ability of the clouds to shield themselves effectively 
against the interstellar radiation field: clouds that are too 
diffuse to shield themselves do not cool and hence form few, 
if any stars. Since effective shielding of the UV background 
is also required in order to form large abundances of H2 or 
CO (Glover et al.|2010| [Glover fc Mac Low||2011| ), this nat- 
urally leads to a correlation between molecular gas and star 
formation without the necessity for a direct causal link be- 
tween them. Our current results do not allow us to establish 
whether the column density of dust required to provide ef- 
fective shielding is independent of the cloud properties, or 
a function of the mean density of the gas; this will be ad- 
dressed in future work. 

Finally, it is interesting to speculate about the possi- 
ble observational consequences of this result. One possible 
way of testing this model would be to look for regions in 
which ongoing star formation is not accompanied by signif- 
icant amounts of molecular gas. If molecular gas is required 
for star formation, then it would be very difficult to explain 
the existence of such regions, whereas they are easily accom- 
modated within our model. Unfortunately, such regions are 
likely to be rare, since they require the molecular gas frac- 
tion to be far out of equilibrium, and as our models show, 
high molecular fractions can be reached within a single free- 
fall given typical Galactic conditions. If our results hold at 
lower metallicities (which is a reasonable assumption, but 
one which must be tested by future simulations), then one 
promising place to look for such regions would be within low 
metallicity star-forming dwarf galaxies such as I Zw 18 or 
DDO 154, since in these systems the characteristic chemi- 
cal timescales will be much longer. Indeed, efforts to date 
to detect CO within these galaxies have been unsuccessful 
( |Leroy et aT1|2007| |Komugi etlu~1|2011| ), although whether 
this is due to a general deficit of molecular gas or simply a 
low CO-to-H2 ratio remains unclear. 
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